Superinflation, quintessence, and nonsingular cosmologies 
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The dynamics of a universe dominated by a self-interacting nonminimally coupled scalar 
field are considered. The structure of the phase space and complete phase portraits are given. 
New dynamical behaviors include superinflation (H > 0), avoidance of big bang singularities 
through classical birth of the universe, and spontaneous entry into and exit from inflation. 
This model is promising for describing quintessence as a nonminimally coupled scalar field. 
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A dynamical system approach to a self-consistent non- 
singular cosmological history is presented in the frame- 
work of the classical Einstein equations with a nonmini- 
mally coupled scalar field. The description of the matter 
content of the cosmos with a single scalar is appropriate 
during important epochs of the history of the universe 
[1]. A crucial ingredient of the physics of scalar fields in 
curved spaces is their nonminimal coupling to the Ricci 
scalar R of spacetime, which is required by first loop cor- 
rections [2] , by specific particle theories [3] , and by scale- 
invariance arguments at the classical level [4]. R is well 
known that nonminimal coupling dictates the success or 
failure of inflationary models [5] ; more generally, it turns 
out to strongly affect the cosmic dynamics, which is qual- 
itatively richer than in the minimally coupled case. We 
show, indeed, that nonminimal coupling leads to new dy- 
namical behaviors, such as a regime that we propose to 
call superinflation (H > 0, which cannot be achieved with 
minimal coupling [6]), spontaneous entry into and exit 
from inflation, with or without a cosmological constant, 
and a possible model for quintessence. Spontaneous su- 
perinflation provides a classical alternative, which is im- 
possible with minimal coupling, to semiclassical [7,8] and 
quantum [9] birth of the universe. 

We consider the nonminimally coupled theory de- 
scribed by the action 



■ g^d^yb -2V + ^Ryb : 



(1) 

where k = 87rG (G being Newton's constant), £ is the 
nonminimal coupling constant, and a cosmological con- 
stant A, if present, is incorporated in the scalar field 
potential V(ip). We use the full conserved scalar field 
stress-energy tensor 



(d a ipd a vj - 2VW)) 



(2) 



(where G M „ is the Einstein tensor), thereby avoiding 

the widespread effective coupling K e s = k (l — k£,ip 2 ) 
in the Einstein equations — kT^. We study 

the dynamics of a spatially flat Friedmann-Robertson- 
Walker universe with line element ds 2 — dr 2 — 
a 2 (r) (dx 2 + dy 2 + dz 2 }. This yields the trace equation 
R = —k(o — 3p), the energy constraint 'AH 2 = na (which 
guarantees that the energy density a > 0), and the Klein- 
Gordon equation, respectively 

6 [1 - £ (1 - 6£) kiP 2 ] (h + 2H 2 ) - k (6£ - 1) vb 1 
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QZuHvjijj - 3H 2 (1 - <V 2 ) +kV = 0, (4) 



dV 
~chp~ 







(5) 



The (time-dependent) equation of state of the vb 
field, rather than being imposed a priori, follows self- 
consistently from the dynamics. The system (3)-(5) re- 
duces to a two-dimensional set of first order equations for 
ip and H = a/a (note that this dimensional reduction is 



not possible for spatially curved universes 
has 

f 



,10]). One 
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where 



3H 2 
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Due to the energy constraint (4), the trajectories are 
restricted to a two-dimensional manifold £ in the 
three-dimensional (^H, ip, ipj phase space, possibly with 

"holes" (dynamically forbidden regions) corresponding to 
Q(H,ip) < (cf. Eq. (6)). S is composed of two sheets 
corresponding to the positive or negative sign in Eq. (6). 
The two sheets smoothly join on the boundary Q = of 
the dynamically forbidden region. 

In the following, for simplicity, we project the dynam- 
ics of the phase space onto the {H, ip) plane, but the true 
nature of £ should always be kept in mind. We now 
restrict to the widely used potential 
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consisting of a mass term, a quartic self-coupling and, 
possibly, a cosmological constant term. For consistency 
with previous works [8] we use the dimensionless symbols 
a = Km 2 /6 (m being the scalar field mass) and to = 
-k 2 A/9. 

The fixed points of the system (3)- (5) are the de Sitter 
solutions with constant scalar field 



Hi 
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6(a — 6£,u>) 
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(n^6a£), and (H,ip) = (±^/-3w/k, o). The fixed 

points (10) exist also for co = 0, due to the pres- 
ence of the matter field t/j (in this case, the two points 

^±^/— 3uj/k, O^j collapse into the Minkowski space fixed 

point (0, 0). We only present here the case of conformal 
coupling, £ = 1/6. 
The function 

LW,i > ) = U 2 + j4> i - — *P 2 + V(i>) (11) 
2 4k 

is such that dL/dt — —3Hip 2 along the trajectories. For 
H > 0, L is a Lyapunov function in a region containing 
the origin; the solutions are then confined by the closed 
lines of constant L, implying asymptotic convergence to 
the fixed points on the H axis. This behavior is con- 
firmed by exhaustive numerical simulations reported in 
the following. We first exclude a cosmological constant 
by setting uj = 0. The phase portrait qualitatively differs 
according to the ratio Q/a. 
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FIG. 1. Qualitative phase portrait for the system (3)-(5). 
Shadowed regions correspond to the dynamically forbidden 
regions (G(H,ip) < 0, cf. Eq. (6)), and the dots to the fixed 
points. Figures a-e, were obtained by using SI/ a, respectively, 
2, 5, 3/2, 3/2, and 1/2, and u> = 0. Figure f corresponds to 
the case u> = —1/10 and SI /(a — uf) = 3/2. 

The case SI = 2a: The Minkowski space (^H,ip,ij?j = 

(0,0,0) is a fixed point, attractive for H > and repul- 
sive for H < 0; the projections of the de Sitter spaces 
(±i?0) ±ipo, 0) are saddle points, i.e. they possess attrac- 
tive and repulsive eigendirections in the phase space as 
shown by the arrows in (Fig. la). The new kind of solu- 
tions, only present in this case Q — 2a and unavoidably 
missed in the approach using k c s(t), 



H(t) 



tanh 
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(where C — H + 2H 2 = —R/6 is constant), corresponds 
to heteroclinic straight lines connecting de Sitter fixed 
points, starting along the repulsive eigendirection of one 
of them and ending along the attractive eigendirection of 
the other (Fig. la). They are tangent to the boundary 
of the forbidden regions at (H,ip) = (0, ±t/>o). For \H\ > 
y/C/2, another straight line solution is obtained from the 
general form 



H{t) = 



2 wie %/^ 



w 2 e 



'C/2i 



,ip = ipo, (13) 



where w\ and u>2 are integration constants. The nonsin- 
gular solutions (12) connect a contracting (t — > — oo) de 
Sitter regime to an expanding one (r — > +oo) and passes 



2 



through a minimum nonvanishing value of the scale fac- 
tor (t = 0). 

In addition to these straight lines, there exist other 
heteroclinic solutions: one starting at (H, ip) = (0, 0) and 
ending at (—Ho,ipo), and another one from (i/oiV'o) to 
(0,0). A third solution starts at the expanding de Sitter 
point and goes to infinity, while another one comes from 
infinity and arrives to the contracting de Sitter point. 
The phase portrait is symmetric about the origin. 

Near the Minkowski fixed point (0,0,0), numerical 
analysis confirms the peculiar behavior suggested by the 
Lyapunov function: orbits approaching this point with 
positive H are attracted to it, bouncing back and forth in- 
finitely many times off the Q = boundary in the (H, ip) 

projection (Fig. la). In the space (H,ip,ip\ these orbits 

are seen to spiral down on a cone towards its apex at the 
origin. The cone results from the union of the two sheets 
in the vicinity of the origin. Along the spiral, the orbit 
passes almost periodically from one sheet to the other 
with period Tbouncc = 27r/m (rbouncc is obtained by an 
asymptotic analysis of the system (3)-(5). Typically, af- 
ter a few bounces, the period coincides with Tbounce with 
good accuracy). A similar behavior for £1 < was re- 
ported in the earlier numerical analysis of [10], but using 
the effective coupling K e ff( T ) and the variables ip and ip. 

In the H < half-plane, the situation is reversed: or- 
bits starting with H < are repelled by the origin and 
depart from it bouncing off the Q = boundary. 

The case £1 > 2a: the situation (Fig. lb) is anal- 
ogous to the previous one, but now the straight hete- 
roclinic solutions are missing, and are replaced by the 
solution starting at the contracting de Sitter fixed point 
and escaping to infinity, and by the solution coming from 
infinity and arriving to the expanding de Sitter point. 
The two-sheeted structure of E implies that no actual 
intersections occur between different orbits in Fig. lb, 
which live in different sheets but are projected in the 
same plane. 

The case a < ft < 2a: as shown in Fig. lc, there 
are no straight heteroclinic lines but new interesting fea- 
tures emerge. A new heteroclinic solution appears start- 
ing from the origin and ending in the expanding dc Sitter 
fixed point. As in the previous case, the quadrant ip > 0, 
H < is obtained from the ip > 0, H > one by reflec- 
tion about the ip-axis and time-reversal. 

The crucial feature of this case is the appearance of 
a dense set of homoclinic solutions (Fig. Id) departing 
from the origin with negative H and returning to it with 
positive H, going around the forbidden region. Super- 
inflation plays a central role along these orbits; only a 
regime with H > permits the smooth transition from 
an initial contracting (H < 0) phase to an expanding 
one (H > 0). This transition occurs at the nonvanish- 
ing minimum of the scale factor. The behavior of this 
family of homoclinics, as well as of the other solutions, is 
universal: they rapidly converge in the spiraling region 
near the origin, irrespective of initial conditions. As all 



of these homoclinics originate around Minkowski fixed 
point due to its instability with respect to perturbations 
with H < 0, this dynamics constitute a classical alterna- 
tive to the previously proposed semiclassical birth of the 
universe from empty space [7,8]. 

The case < < a: the fixed points (10) disappear 
and the only bounded solutions are the homoclinics as- 
sociated with the origin (see Fig. le). This situation is 
therefore the most favorable for the classical spontaneous 
exit from empty Minkowski fixed point. 

Analogous results hold also with a small and positive 
cosmological constant A (see, for instance, Fig. If), with 
the phase portrait being classified according to 0/ (a— ui), 
but the fixed point (0, 0, 0) of the uj = case splits into 
two de Sitter fixed points with memory of the previous 
stability properties. Now, the approximate period be- 
tween two consecutive bounces is 

Tbouncc = , = ■ (14) 
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FIG. 2. The plane (H,ip) for the u = case and the 
equation of state. The darkest region corresponds to the 
superinflation regime (H > 0). In the region H < and 
H > 0, the equation of state associated with the bounc- 
ing solution ip passes periodically through radiation domi- 
nation (crossing the //-axis), matter domination (p — 0), and 
re-acceleration (between Q = and a = lines). Asymptot- 
ically, as t — > +oo, the universe becomes matter dominated 
a(r) oc t 2//3 and tends to infinite dilution. 

Returning to the u> — case, the (H, ifj) plane is di- 
vided into sectors by the straight lines H = ±y/a/2 ip, 
H = ±\/a ip H = ± \J2a ip corresponding, respectively, 
to H = 0, a = 0, and pressure p = (the //-axis 
corresponds to p = er/3). The lines H = ±^a/2ip 
mark the transition between inflationary (a > and 
H < 0) and superinflationary (H > 0) regimes. The 
lines H = ±y/aip divide regions corresponding to infla- 
tion and to decelerated expansion, while H = ±^/2aip 
divide regions of positive and negative pressures (Fig. 2). 
The crucial condition for superinflation to occur is that 
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the line H = (or parts of it) belong to the dynami- 
cally accessible region Q > of the (H, ip) plane. This 
implies that, for an arbitrary potential V(iJ>), superin- 
flation corresponds to ipdV/dip < 0. For our particular 
potential (9), this requires ft > 0. The superinflation- 
ary behavior occurs only once along each homoclinic and 
brings the solution from the primordial Minkowski fixed 
point neighborhood to the succession of eras correspond- 
ing to different equations of state (during each bounce), 
towards infinite dilution and equation of state p = 0. 
Indeed, asymptotic analysis for r — > +oo and for any 
value of a and (with u> = 0) shows that the scale fac- 
tor a(r) exhibits oscillations of concavity corresponding 
to accelerated and decelerated epochs. These oscillations 
are damped as r — > +oo; in this regime, a(r) oc t 2 / 3 and 
the universe becomes matter dominated. 

While it is not claimed here that the evolution of our 
universe is modeled by an entire orbit of the system (3)- 
(5) on the accessible manifold E, the application to spe- 
cific eras of the cosmological history is intriguing. Indeed, 
in the bounces reported above (during which H < 0), 
one encounters, respectively, radiation domination cross- 
ing the iJ-axis, matter domination (p = 0), accelera- 
tion (a possible quintessence model?) until the next 
bounce in the (H,ip) projection, where this sequence is 
reversed. If we identify one period as our cosmologi- 
cal history, then the reported accelerated expansion of 
the universe today [11] suggests to locate our epoch in 
the sector between the line H = ^aip and the Q = 
boundary. The identification of the age of the universe 
(~ I0 17 s) with Tbouncc would then yield the scalar field 
mass m ~ I0~ 13 eV, which is suggestive of an axion 
[1] or of an ultralight pseudo-Goldstone boson (that the 
quintessence field should be very light was already sug- 
gested [12]). 

Although our simplified model is not a quintessential 
scenario, its features are promising for the description 
of quintessence as a nonminimally coupled scalar field. 
Furthermore, the main properties of the model presented 
here remain unaltered for a large class of potentials V(ip). 
Finally, the fact that the dynamics of the system (3)- (5) 
are confined to the two-dimensional smooth manifold E 
indicates the absence of chaos, a conclusion supported 
by further analysis. This is not the case, a priori, of a 
spatially curved universe and/or of multiple matter fields. 
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